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Abstract 

We applied an analytical formula for Bose-Einstein correlations (BEC) devel- 
oped by us recently to high-energy heavy ion collisions, in particular to data on 
S+Pb— > 7r~*'-7r"*"+X reaction at energy 200 GeV/nucleon reported by the NA44 Col- 
laboration. It takes into account both Coulomb and strong (vr-vr s-wave; 1=2) final 
state interactions (FSI). We have found that inclusion of the strong interaction in 
addition to Coulomb correction affects significantly the extracted parameters of the 
BEC like the source size R, the degree of coherence A and the long-range correlation 
parameter 7. In particular, the A parameter of the BEC is increased by about 20%. 
Our results differ from those obtained in e'^e'^ annihilation in the following way: the 
A parameter does not reach 'chaotic limit' and the 7 parameter does not approach 
to zero. 
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Introduction: Recently we have presented an analytical formula for Bose-Einstein cor- 
relations (BEC), which takes into account both Coulomb and strong [tt-tc s-wave; 1=2) 
final state interactions (FSI) and applied it to e^e~ annihilation data [jl|. Our formula has 
been constructed by using a wave function obtained from applying the phase shift of the 
strong interaction to the asymptotic Coulomb wave function, and extrapolating it to the 
non-asymptotic region with the normalization of the Coulomb wave function. We shall 
apply it here to analysis of high energy nuclear collisions using in addition the following 
important parameters entering the BEC: degree of coherence A, long range correlation 7 
and normalization c. This will make our analysis of the FSI much more clear when com- 
pared with a 'standard' fitting formula with or without the Coulomb corrections. In the 
following sections we first provide brief explanation of our formula, use it next to inves- 
tigate data for S+Pb— 7r+-7r+ +X collisions at energy 200GeV/nucleon reported by the 
NA44 Collaboration H] and, finally, we close our discussion with some concluding remarks. 



Coulomb and strong wave function with s-wave phase shift: Let us consider the 
Coulomb asymptotic wave function for the identical tt-tt scattering (of relative momentum 
Q, the pion mass is m^r). Its form is well known and given by 0, 

n2 



\E'^^-^™(k, r) = exp{i{kz + riln{k{r — z)))} 1 + 
, /■//)Aexp{«(fc?^ -r]\n{2kr))} 



ik{r — z) 



where Q = 2k, z = rcos9, r] = ma/Q and the scattering amplitude f{9) is given by: 

rj 1 



m 



exp{— 2z?7 lnsin(^/2) -|- 2i argr(l + if])}. 
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The s-wave component of the \E^Q'*'™(k, r) wave function is given by the following for- 
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+ ex'p{—i{kr — ?7ln(2A;r))} — 
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The strong wave fucion of the s-wave has been obtained by extracting the out-going 
wave (up to the order 0{l/r)) of the eq.(^ and using the phase shift of the tt-tt scattering 



induced by the strong interaction 0-|TO|. This phase shift is phenomenologically given 
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by (g: 

''^'-1(tT^) (-1.5 < < -0.7 (GeV-) ). (3) 

As a result we obtain the following strong interaction wave function in the asymptotic 
region: 

^..(k.r)^/°W ''''P'''*^'--"'°P*^'''» . (4) 

r 

fi0) = ^exp(2zr/o)(exp(2z5f ) - 1). 

Therefore, the total asymptotic wave function is expressed as the following sum of both 
components: 

v^totai(k, r) = vl/-y-(k, r) + r). (5) 

In the small kr (internal) region the Coulomb wave function is given by Q 

*c(k,r) = T{l+i7])e-'''^/^e^-''F{-tr];l;ikr{l-cose)), (6) 

(where F denotes the confluent hypergeometric function) whereas the exact form of the 
strong wave function in this region is unknown. We assume therefore that the strong wave 
function is given by the extrapolation of the asymptotic wave function into this internal 



region |12] with the normalization of the Coulomb wave function |jl|]: 



vi>totai(k, r) = xi>c(k, r) + ^G{r])MK r), (7) 

where G{ri) denotes the usual Gamow factor: G{r]) = 27r?7/(exp(27r?7) — 1) (for the discus- 
sion on the validity and usefulness of this assumption and on the smooth connection of 
wave functions given by eqs.(^ and (|^ cf. Ref. ||l|]). Because in the numerical calculations 
of the ^c(k, r) one encounters a wild oscillating behavior in larger kr region, we use here 



the 'seamless' fitting method developed in [p 



Theoretical formula of BEC with Coulomb and strong wave function: To 

describe a pair of the identical bosons, we have to symmetrize its total wave function in 
the following way: 

Au = ^[^c(k, r) + vl/g(k, r) + <|.,t(k, r) + <|.f,(k, r)], (8) 

where the superscript S denotes the symmetrization of the corresponding wave function. 
Here $st(k, r) stands for the wave function induced by the strong interactions. Assuming 
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that a source function is given by p(r) we obtain the following expression for the BEC 
including the FSI pi]: 



(9) 



oo 
m,n=0 



m + n + 1 



Iri{2 + m + n)Ai{n)A\{m) x 



1 + 



n\m\ 



[n + m)\ 



n 



17] 



1 - 



17] 



= (1 + Aic) + (^2B + Aec 

Icst = 23? 



o oo 

{2ky^ exp {-i{r]o + 5^ )) sin 5^'^ ^ ^R2(l + ri)A2(n, 0) 
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-E'2B 

1 + Aic 



d^r p{r) e 



iQ r 



^0 (2n + l)!(2n + l)(2n + 2)^ 



T{7i + n]) {-2ikY 
T{i7]) (n!)^ ' 

r(n + / + 1 + i7])T{2l + 2) (-2iA;)" 
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47r J drr'^p{r), 



A,{7l) 

A2{n,l) 

In this paper N^"^"^^ /N^^ stands for the FSI corrected ratio of pairs of the identical 
charged bosons measured in a single event to a product of single bosons taken from 
different events. The Aj^q and AgQ originate from the correction for the finite size effects 
in the Coulomb interaction. In the following we shall use an artificial Gaussian form of 



the source function: p(r) 
Fourier transform: 



1 \3 



exp 



. 2/32 / ' 



which results in the following form for its 



^2B = exp(-/5'Q72). 



Reanalyses of data reported by NA44 Collaboration: The parameters mentioned 
at the beginning, namely: degree of coherence A, long range correlation 7 and normaliza- 
tion c, are introduced into our formula eq.® in the usual way leading to the following 



final formula: 



^(±±)/^BG^g^2A;) = C (1 + Alc + Aec + /est + /st) 

E2B 



X 



1 + A- 



:i + 7Q). (10) 



1 + Aic + Aec + /est + /st- 
This formula have been applied to data on S+Pb vr+vr"'" +X reaction at energy 200 
GeV/nucleon reported by the NA44 Collaboration 0. 

Before presenting our results based on eq.(|10|) we must, however, first clarify one 
point concerning aplication of the Coulomb correction. We show the results in Table 
I obtained by applying the following pure Coulomb (i.e., without strong interactions) 
corrected formula: 

E2B 



^(±±)/^BG = c (1 + Aic + Aec) x 



1 + A- 



:i+7Q), (11) 



1 + Aic + Aec- 

to the BEC data presented by the NA44 Collaboration. On the other hand, the NA44 
Collaboration have also corrected their data by Coulomb correction factors, K^^^i'^, which 
have been calculated by using the Coulomb wave function integration. However, as is 
shown in the Table I, slight discrepancies are found between their results (obtained by 
using the K^^^^) and ours (obtained from minimum fit using formula eg. ([TT|) ) Q The 
origin of discrepancy shown in the Table I can be traced down to the 15% differences 
between the K^^^^ and the K°^J^i (see Table II) (one can see this by evaluating the 
correction factor, K^^^^ = 1/G{ri)[l + Aic + Aec])- Out final results obtained by using 
eq.(p!OD are shown in Tables III, IV and in Fig.l. For the sake of reference the results of 
eq.(pJ]) and of the 'standard' fitting formula: 

are also shown there. It is found that the FSI corrections, which take into account not 
only the Coulomb but also the strong interactions between the identical charged pions, 
affect substantially the source size parameter R, the degree of coherence parameter A and 
the long-range correlation parameter 7. When compared with the Coulomb correction 
case only (i.e., with eq.(ll)), the R parameter is reduced by about 5% and the A param- 
eter is increased by about 20% remaining, however, always below its 'chaotic' limit value 



^To compare our estimated values with those reported by the the NA44 Cohaboration, we have to 
use a relation, /3 — \plR ( which reflects the fact that NA44 Collaboration employed the following fitting 
formula: c (1 + A exp(-g2i?2)) ). 
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of A = 1[||] (notice that in e+e~ annihilation the corrections caused by the FSI pushes the 
A towards this hmit). In what concerns the parameter 7, we have found that it does not 
reduce to zero (differently than in the e^e~ annihilation case). 

Concluding remarks: We applied our improved formula for the BEC (which accounts 
for the FSI) to data for S+Pb 7r+7r+ +X at energy 200GeV/nucleon presented recently 
by the NA44 Collaboration. We have found that inclusion of the full FSI (including 
also the strong interactions, cf eq.(|TUp) increases the extracted value of the degree of 
coherence parameter A by about 20% in comparison to the case where only the Coulomb 
correction are applied (cf. eq.(0)). However, differently than in e^e~ annihilation case 
(cf. [|l|) it does not reach 'chaotic' limit of A = 1. We have also found that the long-range 
correlation parameter 7 does not approach to zero which also seems different from the 
e^e~ annihilation case. On the other hand the source size parameter R is decreased by 
about 5% only. 
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Figure Captions. 

Fig. 1. Analysis of data for S+Pb^ vr+vr"'" +X reaction at energy 200 GeV/nucIeon re- 
ported by the NA44 Collaboration [Q] by the minimum fitting method. The strong 
interaction phase shift parameter has been set to oq = —1.20 GeV^^. The sohd and 
dotted curves represent fitted result obtained by using eq. ([TO|) with and without 
factor (H-7Q), respectively. 

Table Captions. 

Table 1. Results of analysis of BEG data for S+Pb^ 7r+7r+ +X reaction at energy 200 
GeV/nucleon reported by the NA44 Collaboration [0]. Comparison of our results 
obtained from fitting of the eq.(ll) with values reported by the NA44 Collaboration 
which were obtained by using the Coulomb correction factor /T^^^/IQ]. 

Table 2. Comparison of the Coulomb correction factor K^^^^ (used in ref.[2]) with 
the K^J^i (evaluated from eq.(2)). 

Table 3. Results of analysis of BEG data for S+Pb^ 7r+7r+ +X reaction at energy 200 
GeV/nucleon reported by the NA44 Collaboration [^]. Comparison of our results 
obtained from fitting of eqs.(lO), (11) and (12) without factor (1 + 7Q). We have 
employed the parameter ao = —1.20 GeV~^. 

Table 4. Results of analysis of BEG data for S+Pb-^ tt+tt^ +X reaction at energy 200 
GeV/nucleon reported by the NA44 Collaboration [0]. Comparison of our results 
obtained from fitting of the eqs.(lO), (11) and (12). We have employed the parameter 
ao = -1.20 GeV-^ 
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Table I: 



formula 


c 


R [fm] 


A 


7 


XVNDF 


correction by 

7>'NA44 
^Coul 


1.000±0.002 


4.50±0.431 


0.46±0.04 




18.1/16 


cq. (11) 
without (1 + 7(5) 


1.000±0.003 


4.697±0.272 


0.481±0.031 




34.5/16 


eq. (11) 


1.031±0.009 


5.159±0.331 


0.453±0.034 


-0.181±0.048 


22.1/15 



Table II: 



[MeV/c] 


^NA44 
Coul 


Coul 


5 


1.565 


1.807 


15 


1.161 


1.112 


25 


1.074 


1.037 


35 


1.037 


1.014 


45 


1.020 


1.008 


55 


1.012 


1.006 


65 


1.009 


1.006 


75 


1.007 


1.006 


85 


1.005 


1.005 


95 


1.005 


1.004 


110 


1.003 


1.003 


130 


1.003 


1.002 


150 


1.002 


1.000 


170 


1.002 


0.999 


190 


1.001 


0.999 


225 


1.001 


1.000 


275 


1.001 


1.003 


325 


1.001 


1.004 


375 


1.000 


1.007 
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Table III: 



formula 


c 


R [fm] 


A 


7 xVndf 


eq. (10) 

strong + Coulomb 


1.007±0.003 


4.422±0.282 


0.593±0.028 


— 46.5/16 


eq. (11) 

Coulomb 


1.000±0.003 


4.697±0.272 


0.481±0.031 


— 34.5/16 


eq. (12) 

standard 


1.019±0.003 


4.209±0.206 


0.528±0.031 


— 48.1/16 



Table IV: 



formula 


c 


R [fm] 


A 


7 


xVndf 


eq. (10) 

strong + Coulomb 


1.045±0.008 


4.870±0.320 


0.541±0.031 


-0.232±0.047 


24.6/15 


eq. (11) 

Coulomb 


1.031±0.009 


5.169±0.331 


0.453±0.034 


-0.181±0.048 


22.1/15 


eq. (12) 

standard 


1.069±0.010 


4.929±0.292 


0.514±0.034 


-0.274±0.051 


18.1/15 
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1 1 1 


1 1 1 1 

eq.(10) without (1+yQ) 


o 




eq.(10) 






^ NA44 Collab. 
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